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Abstract 

Several methods are currently available to simulate paths of the Brow- 
nian motion. In particular, paths of the BM can be simulated using the 
properties of the increments of the process like in the Euler scheme, or 
as the limit of a random walk or via L 2 decomposition like the Kac- 
Siegert/Karnounen-Loeve series. In this paper we first propose a IFSM 
(Iterated Function Systems with Maps) operator whose fixed point is the 
trajectory of the BM. We then use this representation of the process to 
simulate its trajectories. The resulting simulated trajectories are self- 
affinc, continuous and fractal by construction. This fact produces more 
realistic trajectories than other schemes in the sense that their geometry 
is closer to the one of the true BM's trajectories. The IFSM trajectory of 
the BM can then be used to generate more realistic solutions of stochastic 
differential equations. 
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1 Introduction 

In this paper we show how to solve the inverse problem for IFSM in the case 
of trajectories of stochastic processes in L 2 . The method is based on the solu- 
tion of the inverse problem for IFSM due to Forte and Vrscay [3]. This is an 
extension of classical IFS methods which can be used for approximating a given 
element of L 2 (H) thus in particular trajectories of stochastic processes on this 
space. The final goal of this approach is simulation. Indeed, several methods are 
currently available to simulate paths of stochastic processes and in particular of 
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the Brownian motion. Paths of the BM can be simulated using the properties 
of the increments of the process like in the Euler scheme [B] , or as the limit of 
a random walk or via L 2 decomposition like the Kac-Siegert/Karnounen-Loeve 
series |7]- In this paper we first propose a IFSM (Iterated Function Systems 
with Maps) operator whose fixed point is the trajectory of the BM. We then 
use this representation of the process to simulate its trajectories. The result- 
ing simulated trajectories are self-affine, continuous and fractal by construction. 
This fact produces more realistic trajectories than other schemes in the sense 
that their geometry is closer to the one of the true BM's trajectories. The IFSM 
trajectory of the BM can then be used to generate more realistic solutions of 
stochastic differential equations. The paper is organized as follows: Section 
recalls the theory of IFSM on L 2 , Section [3] recalls some details for stochastic 
processes with trajectories in L 2 and the link with the IFSM theory. Section 0] 
presents the application of the IFSM theory to the problem of simulation with 
particular attention to the case of the Brownian motion. 

2 IFS with Maps (IFSM) on L 2 {H) 

The basic idea of Iterated Function Systems (IFS) can be traced back to some 
historical papers but the use of such systems to construct fractals and other 
similar sets was first described by Hutchinson (1981). The fundamental result 
on which the IFS method is based is Banach theorem. The mathematical context 
is the following: given y in a complete metric space ( Y, d) , find a contractive 
operator T : Y — > Y that admits a unique fixed point y* G Y such that dy(y, y*) 
is small enough. In fact if one is able to solve the inverse problem with arbitrary 
precision, it is possible to identify y with the operator T which has it as fixed 
point. The fundamental theorems on which the IFS method is based on are the 
following: 

Theorem 2.1. (Banach Theorem) Let (Y, dy) be a complete metric space; sup- 
pose there exists a mapping T : Y — > Y such that 

d Y (T(x),T(y) < cd Y (x,y) 

for all x,y £ Y and some c € [0, 1). c is said to be the contractivity factor of 
T . Then there exits a unique y* 6 Y such that T{y*) — y* and for any y £ Y 
we have d(T n (y),y*) — > when n — > +oo. 

Theorem 2.2. (Collage Theorem) Let (Y,d Y ) be a complete metric space. 
Given y € Y suppose that there exists a contractive map T with contractiv- 
ity factor c G [0, 1) such that dy{y 1 T{y)) < e. If y* is the fixed point of T then 
dy(y,y*) < ^- c . 

Theorem 2.3. Let (Y, dy) be a complete metric space and T%,T2 be two con- 
tractive mappings with fixed points y\ and y% ■ Then 

d Y{vt,yl) < -7— — d YtSup (T 1 ,T 2 ) 
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where 

dY,sup{Ti,T 2 ) = supd(T 1 (x),T 2 (x)) 

x<£Y 

and C\ is the contractivity factor of T\ . 

We are going to use a particular class of IFS operators, known as IFSM (IFS 
with Maps), introduced by Forte and Vrscay in 1994. Let H = [0, 1], /z be the 
Lebesgue measure on B(H) (the Borel cr-algebra) and for any integer p > 1 
let L P (H) denote the linear space of all real valued functions u such that u p 
is integrable on (B(H),^). To build a contraction map T on L 2 (H) we need 
TV-map contractive IFS i.e. a set of maps W = {w\, w 2 , ■ ■ • , w n } and a set of 
functions (grey level maps) <j> = {<p%, cj) 2 , . . ■ , 4> n } with : R — » M. The operator 
T corresponding to the N map IFSM(ui,(/>) is 

(T«)(a:) = X;>*("(«'* C 1 ) 
fc=i 

where the prime means that the sum operates on all those terms for which w^ 1 
(the inverse function of Wk) is defined. Let us define the following two sets 

Sim(H) = {w:H^H:3ce [0, l),\w(x) - w(y)\ = c\x-y\,Vx,y<E H} 

£ip{R) = {<f>: R^ M : 3K G [0, oo), \</>(ti) - <f>{t 2 )\ < K\t x - t 2 \,Mt u t 2 G M} 

Theorem 2.4. Let (w, 4>) be an IFSM such that Wk G Sim(H) and <j) k G 
££p(R) /or 1 < K < N. Then T : L 2 {H) -> L 2 (#) and /or any u,w G L 2 (H) 
we have 

d 2 (Tu,Tv) < Cd 2 (u,v) 

where 

N 

fe=i 

Given u G L 2 (H) the inverse problem consists of finding the operator T such 
that 

N 

u(x) = (Tu)(x) = ^2'(j)k{u(w^ 1 (x))). 

k=l 

In |3] it is proved that this problem can be reduced to the determination of grey 
level maps 4>k which minimize the collage distance A 2 

A 2 = \\v~Tv\\ = 

N 

/ ^'HkH'Wk'i^-^W- 
Jh k =i 

In the special case when 
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• U*=i H k = Ufc=i Wi{H) = H i.e. the sets H k "tile" H 

• n(wi{H) n Wj(H)) = for i ^ j 

we say that the maps u>fc are nonoverlapping. Later, in the applications, we will 
assume that, for 1 < k < N, 

• w k (x) = s k x + a k 

• < c fe = \s k \ < 1 

• <p k (t) = a fc i + k , K k = \a k \ 
The collage distance becomes 

A 2 =< t> — Tu, v — Tv >= 

AT N 

< ip k ,ipi > a k ai + 2 < ip k ,ipi > a k (3i+ < £ fe ,& > 

N 

-2 < t> , ipk > a k + <v,£, k > f3 k + < v, v > 
fc=i 

where 

^fe(x) = u(tu^ 1 (a;)), £ k (x) = I Wk ( H )(x) 
A 2 is a quadratic form in aj and /3j, that is 

A 2 = .t t Ar + 6 T a: + c (2) 

where x = (ot\, . . . a k , . . . , (3 k ). The matrix A is symmetric and 

a i,j =< tyutyj >,0>N+i,N+j =< £i>£j > 

a i>N+ j =< ipi,£j >,bi = -2 < u,^ >,&Ar +i = -2 < > 

and c = As in [5j we add an additional constraint in order to guarantee 

that the minimum of this quadratic form exists on a compact subset of feasible 
parameters on and The additional constraint is 

N 

IMIi + /3fc) -Nk<o. 

fe=i 

The maps Wfc are choosen in an infinite set W of fixed affine contraction maps 
on iJ which has the /i-dense and nonoverlapping property (in the sense of the 
following definition); When (a k , (3 k ) — (0, 0) the corresponding w k is superfluous 
and the fc-th term can be dropped from Q). 

Definition 2.1. We say that W generates a \x-dense and nonoverlapping family 
T of subsets of H if for every e > and every B C H there exists a finite set 
of integers i k , i k > 1, 1 < k < N , such that 
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• A = U% =1 w ik (H) C B 

• (jl(B\A) < e 

• l i(w ik (H)nw il (H)) = Oifk^l 



Let 



W N ={ Wl ,...w N } 



be the N truncations of w. Let § N = {0xj • • • j <Aiv} the TV vector of affine 
grey level maps. Let x N be the solution of the previous quadratic optimization 
problem J2J) and Ajy min = A 2 N (x N ). It can be shown that min may be 
arbitrarly small when N — » oo (see [3]). 



Let (fi, JF, P) be a probability space and {J-"t, t G H} be a sequence of cr-algebras 
such that Tt C T. Let : fl x i? — > R be a stochastic process in L 2 (H). 

that is a sequence of random variables J^-adapted (that is each variable X(tu,t) 
is J^-measurable). Given u> G a trajectory of the process is the function 
X(w,t) : H — > R belonging to L 2 (H). For a given the trajectory of 

the stochastic process, the aim of the inverse problem consists in finding the 
parameters of the IFSM such that X(u>,t) is the solution of the equation 



In this case the coefficients of the matrix A and the vector b of the previous 
section become 



3 IFSM for stochastic processes on L 2 {H) 



X(oj, t) = TX(u, t) for a.a. uj e fl 



chj {uj) 



L 




and if i = j it becomes 




The other elements in the matrix A can be calculated as 




and 




For the vector b 



bi{u>) 



2<X,ipi> 



L 



X{u,t)X{w,wr x {t))dt 
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and 

b N+l (cj) = -2<X, & >= f X(lj, t)dt 

JwiiH) 

In the nonoverlapping case, we have 

• dij = 0, i 7^ j, and aj ; j(u;) = Cj J H X 2 (u; 7 t)dt, 1 < i,j < N 

• a N+i , N+ j = 0, 1 < i, j < TV, i ^ j and a N+itN+i = fi(wi(H)) 

• ai,Af+j = 0, 1 < i, j < N, i ^ j and a itN+i (uj) = Ci f H X(w,t)dt 
It also holds this self-similarity property. 

Theorem 3.1. Let (ak,f3k) be the solution of the inverse problem with a set of 
nonverlapping maps u>k and suppose that X(u),t) = TX(u,t). Then 

X(u, u>i(t + h)) - X(w, u>i(t)) = ai(X(u}, t + h)- X(oj, t)). 

for alll<i<N. 

Proof. In fact we have 

X(oj, w t {t + h)) - X(oj, Wi{t)) = TX(oj, Wi(t + h)) - TX(u, w,(t)) 

N N 

= a k (X{cj, w^ 1 (w l {t + h) j) +f3 k -J2 a k(X{tJ, w^ 1 {w l {t) j) + [3 k 

k=l k=l 

= a t {X{ij 1 wl 1 {w l {t + h))) + <3i- aiiXiu^^iwitt))) + fa 
= a,(X(uj,t + h) -X(uj,t)). 

□ 

3.1 The Kac-Siegert decomposition of L 2 (H) stochastic 
processes 

We suppose that a.e. X(u,t) is an element of a subspace S of L 2 (H) and that 
X(ui, t) is a zero-mean process. Let K be the covariance function of this process 
that is 

K (s, t) = Cov [X(u, s), X(u, t)]. 

and assume 

/ K(t, t)dt < oo. 
iff 

If Ai > A 2 > ... > comprises the entire spectrum of eigenvalues of K, where 

/ f(s)K(s, t)dt = Xf(t), 0<t<l 
Jh 
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and the associated ortho-normal cigcnfunctions fi form a complete set of the 
subspace S then the Kac and Siegert decomposition holds: 

oo 

K(s,t)=Y,\ j f j (8)f j (t), 0<s,t<l 
i=i 

We also have that 

Zi = f X(t)fj(t)dt 

J H 

are uncorrelated random variables with mean and variance . The following 
theorem states some properties of this decomposition. 

Theorem 3.2 (see Ch.5, |7])- Suppose that X and K satisfies the properties 
above. Then 

• i) E£L X i < 00 

• Ej=i Zjfj ~*\\\\2 X as m — > 00 a.s. 

• Hi) Zj =< X, fj > are with mean and variance \j 

. iv) J H X*(t)dt = J2?=i Z] = ET=i W ™here Z* = *L- 

• v) E[X(t) — Y^jLi Zjfj(t)] 2 — > for each t as m — > 00 

• vi) X = E^Li \/~\ifjZj with Z* = —J= uncorrelated with mean and 
variance Xj . 

4 Simulation of Brownian motion via IFSM on 

= [0, 1] 

In the literature there several methods of simulation of the trajectory of the 
Browian motion, i.e. the stochastic process {B(uj,t),t G [0,1]}, such that 
-B(O) = a.s., B(t) — B(s) is distributed with Gaussian law with zero mean 
and variance t — s, and with independent increments. 

The Euler method 

In this case, the trajectory is obtained simulating the increments of B in the 
following way: B(0) = 0, B(t i+1 ) = B(ti) + \/U+i — U ■ Zi, where the ZiS are 
independent 7V(0, 1) random variables. In the other points the trajectory is 
built by linear interpolation of these simulated data. 
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The Kac-Siegert method 



Karhunen-Loeve / Kac-Siegert decomposition of B is better for pathwise simu- 
lation 



4>i a basis of orthogonal functions and Z^s are jV(0, 1) 

The trajectory generated by Euler method is too simple and regular to mimic 
the roughness of the BM; moreover the simulated path is stochastically equiva- 
lent to the true trajectory only on the points of the grid used in the simulation. 
The Kac-Siegert decomposition of the BM is a pathwise approximation which 
can lead to a too smooth path (see figure QJ; our idea is to use IFSM for gen- 
erating fractal trajectories of the BM. There are applications in finance (for 
instance pricing of american options) in which the whole path matters; our 
IFSM approch produces a global approximation of the trajectory preserving the 
geometric fractal nature of the target. 

This method can be also used to simulate paths of solutions of stochastic 
differential equations driven by Brownian motion (e.g. diffusion processes) re- 
placing the linear behaviour of the Euler trajectory with a fractal object. For 
the solution of the inverse problem for the BM, we choose the so-called wavelet 
type maps jo], that is: 



with i = 1,2,... and j — 1 ... 2* For each fixed i, the family w*j is a set of 
nonoverlapping maps. For these maps = 2~ l < 1. We organize them as 
follows 



To simulate a trajectory of B with non overlapping maps we then need to 
simulate the joint distribution of all this objects 



B(u,t) =^Z i <j> i (t), 0<i<l 



with 




W\ = W^i W2 — U>* 2 W3 — W21 W4 — W* 2 2 ■ ■ ■ 





I 



1 



B{t)B{wr l {t))dt= [ B(t)B 
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dt 



Jo 
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Figure 1: Paths of Brownian motion simulated by the Euler scheme (dotted 
line) and using Kac-Siegert decomposition (continuous line). The same (n = 25) 
pseudo-random Gaussian numbers were used. 



but it appears to be still a too difficult problem. 

In practice, it is preferable to use all the above maps and not only the subset 
of non-overlapping maps. In this case, we need simulate the value of the trajec- 
tory of the Brownian motion on a fixed grid (using one of the known methods) 
and we use these points to approximate the integrals in the quadratic form. We 
then solve the constrained quadratic programming problem using standard algo- 
rithms (see e.g. 0). Figure [21 bottom shows an example of trajectory generated 
using the IFSM approach using wavelet type maps for i = 1, . . . , M, M = 8. 
Figure 0top represents the Euler trajectory built on 50 Gaussian terms which 
has been used to build the IFSM. As one can notice the IFSM path shows more 
"fractal" complexity then the corresponding Euler path. 

5 Conclusions 

We have proposed a new method to generate paths of the Brownian motion. 
These IFSM paths seems to mimic more closely the fractal nature of the tra- 
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jectory of the Brownian motion then existing schemes. At current stage we 
are not able to show formal property of the IFSM path in terms of strong 
and weak approximation (see [B]). Open source software for generating IFSM 
trajectories written in C and R language [S] is available via ifs package at 
http : //CRM . R-pro j ect . org for free download. 
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